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Abstract. We report an unexpected systematic degeneracy between different 
multiplcts in an inversion symmetric system of two coupled Gaudin models with 
homogeneous couplings, as occurring for example in the context of solid state 
quantum information processing. We construct the full degenerate subspacc 
(being of macroscopic dimension), which turns out to lie in the kernel of the 
commutator between the two Gaudin models and the coupling term. Finally we 
investigate to what extend the degeneracy is related to the inversion symmetry 
of the system and find that indeed there is a large class of systems showing the 
same type of degeneracy. 
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1. Introduction 

In a large variety of nanostructures spins couple to a bath of other spin degrees of 
freedom. Commonly such systems are described by so-called central spin models. 
Important examples are given by semiconductor [TJ [21 121 H] and carbon nanotube [SJ 
quantum dots, phosphorus donors in silicon [6J, nitrogen vacancy centers in diamond 
|2l H I] and molecular magnets [10]. Motivated by the perspective to utilize the 
central spins as qubits [TTJ [12] or to effectively access the bath spins via the central 
spins [131 [H], presently central spin models are the subject of extensive theoretical 
as well as experimental research. However, their importance in a more mathematical 
context became clear already several decades ago, when Gaudin proved the Bethe 
ansatz integrability of the central spin model with one central spin (Gaudin model) 
[15]. Since then they are in the focus of the important field of quantum integrability 



It is well-known that the energy levels of a quantum system usually tend to repel 
each other and degeneracies are exceptional events [22^ . Hence there are only extremely 
few examples of systems with degenerate eigenstates and even less, whose eigenstates 
are systematically degenerate. Famous examples are given by the hydrogen atom [23 , 
the n-dimensional harmonic oscillator [24] or the Haldane-Shastry model [25] [26] . In 
all three cases the degeneracies are due to hidden symmetries requiring a dedicated 
analysis. 

In the present work we we study the spectrum of an inversion symmetric central 
spin model consisting of two coupled Gaudin models with homogeneous coupling 
constants, meaning they are chosen to be equal to each other. In order to lower 
the dimension of the problem, the baths of the two Gaudin models are approximated 
by single long spins, which does not change the set of eigenvalues of the Hamiltonian. 
Surprisingly, the resulting Hamiltonian exhibits systematically degenerate multiplcts 
of consecutive total angular momentum and alternating parity, a situation somewhat 
similar to the degenerate multiplets of orbital angular momentum in the hydrogen 
atom. 

The outline of the paper is a follows. The degeneracies in the coupled Gaudin 
models are first analyzed in a numerical approach in Sec. [2] In Sec. [3] we analytically 
construct the full subspace of degenerate states which turns out to be located in the 
kernel of the commutator between the two Gaudin models and the coupling term. 
In Sec. [4] we furthermore investigate the role of the inversion symmetry and show 
that indeed there is a whole class of systems with spectra showing the same type of 
degeneracy. 

2. Model and spectral properties 

The Gaudin model [TS] describes the coupling of a central spin Si to a set of bath 



via some coupling constants A\ , which have the unit of energy. In the following we 
choose the couplings to be homogeneous, i.e. A\ = Aj. In this case the central spin 



[m nil ng ng eqi ei] . 



spins I\ 
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couples to a simple sum of spins, denoted by Ii from now on. Furthermore we assume 
Si = 1/2. Coupling together two such Gaudin models Hq := Hi + H 2 by 

H c = JexS\S 2 (2) 

finally yields the Hamiltonian subject to our investigation: 

H = H G + H C = A x §ih + A 2 S 2 I 2 + Jex5i5 2 . (3) 

The parameter J ex can be viewed as an exchange coupling. Obviously the Hamiltonian 
conserves the total spin J = S\ + S 2 + I\ + I 2 as well as 7f and i| . 

The ni bath spins couple to different values Ii of the total bath spin squared. 
In the following we study the spectrum of the Hamiltonian for A\ = A 2 = A, where 
A := (1/2)(j4i + A 2 ), on subspaces I\ = I 2 =: I. On these subspaces, in addition 
to the symmetries mentioned above, H is invariant under "inversions" , meaning an 
interchange 1h2. It is clear that this is not the case globally, i.e. on the entire Hilbert 
space. However, subspaces with Ii = I 2 lie fully in the kernel of the commutator [H, r], 
where r denotes the inversion operator. Obviously this only has the two eigenvalues 
(±1). In the following we refer to this as positive and negative parity. 

In order to reduce the dimension of the problem, we approximate each bath by one 
single spin of length /. This neglects the quantum numbers associated with a certain 
Clebsch-Gordan decomposition of the respective bath and therefore the multiplicity 
of the eigenvalues changes, but not the set of eigenvalues itself. Every energy in the 
resulting spectrum indeed appears x±x 2 times in the spectrum of H, where Xi denotes 
the number of multiplets with the quantum number Ii. If for example I\ = i, we 
have E3 



m 

- Id V^-h-l 



(4) 

2 J V \ 2 ll V - 

However, it should be stressed again that the energy eigenvalues themselves remain 
unaltered. 

In Figs. [T] and [2] we show spectra obtained numerically for different values of the 
exchange coupling constant J ex , both for an even and an odd /. Although the spectra 
are quite rich in detail, their global structure becomes already plausible from simple 
qualitative arguments. Obviously we always have four "branches" of energy levels, 
where, in particular for large J ex , three of them form a bundle separated from the 
fourth one. The three former branches consist of states where the two central spins 
are predominantly coupled to a triplet (which has eigenvalue J ea; /4 under H c ) while in 
the latter branch the central spins are mainly in the singlet state (having eigenvalue 
— 3J ex /4 under H c ). The coupling of the central spin triplet and singlet to the bath 
spins leads then leads to the observed further energy splittings between and within 
the corresponding branches. 

An unexpected particular feature, however, occurs in the triplet branch of 
intermediate energy. Here all multiplets are energetically completely degenerate with 
eigenvalue (J ex — 2A)/A. These multiplets have consecutive total spin between J = 1 
and J = 21 and alternating parity. Here positive (negative) parity corresponds to 
21 — J being even (odd). The latter observation is reminiscent to the degenerate 
multiplets of orbital angular momentum found in the hydrogen problem. In general, 
such systematic degeneracies are extremely rare, and hence our finding is interesting 
on its own right. Moreover, since the degenerate subspace is of particularly high 
dimension, potential applications in, for example, solid state quantum information 
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Figure 1. Spectrum of H with baths approximated by two single spins of equal 
length / = 20 for different values of the exchange coupling J ex . The energies (in 
units of A) are plotted against the total spin J, i.e. each data point represents a 
multiplet of 2 J + 1 states. States of positive (negative) parity are signalled by a 
cross (circle). For all exchange couplings Jex 

we have four "branches" of energy 
levels, where the above three ones originate from triplet states with respect to 
H c and the lower one is associated with the singlet state. The triplet branch 
of intermediate energy consists completely degenerate multiplets of alternating 
parity. 
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Figure 2. Analogous data as in Fig. [l] for an odd spin length I = 21. Again we 
find a completely degenerate triplet branch. 
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processing can be envisaged: It is clear that states with overlap exclusively in a 
degenerate subspace do not show any non-trivial time evolution. Therefore, such 
spaces have the potential to provide valuable implementations of long lived quantum 
memory, where the present one appears to be particularly suitable due to its enormous 
size. Note that even in the "thermodynamic limit" J> la fourth of the Hilbert space 
is degenerate: The dimension of the full Hilbertspace is 4(27 + l) 2 and the degenerate 
subspace Hd has dimension J2 n =i(^ n + 1) = + 1)' yielding 

1(1 + 1) _ 1 

(2/+1) 2 ~ 4' 

if I 3> 1. Furthermore, the space of degenerate states detected here decomposes 
into subspaces of different parity which could also serve as a computational basis for 
quantum information processing. 



3. Construction of the degenerate subspace 

So far we have reported on numerical observations revealing an unexpected systematic 
degeneracy in the spectrum. In the following we analytically construct the subspace 
Hd of these degenerate multiplets. 



3.1. General ansatz and first consequences 

As we shall see below, the degenerate states are simultaneous eigenstates of the Gaudin 
part Hq of the Hamiltonian and the coupling H c between the central spins. In other 
words, %d lies entirely in the kernel of the commutator 



[H G , H c ] = -iAJ ex (Si x S 2 j • (Jr - I 2 j . (5) 

Let us first turn to a single Gaudin Hamiltonian, Hi = AS-Ji, with Si = 1/2 and 
Ii = I. The eigenvalues read 

E± (A,I) = ^(±(l + l)+l) , (6) 



2 V V 2 / 2, 

and the eigenstates are given by a well-known Clebsch-Gordan decomposition [55] 

\I±\,m l ) = /^(mOI W,m t - 1) ±M*K)| i)\I,rrn + 1), (7) 
where, apart from standard notation, we have introduced 



The eigenvalues of Hq = Hi + H 2 now follow immediately 

H G \+,mi)\+,m 2 ) = -{h + h)\+,m 1 )\+,m 2 ) (9a) 

# G |+,mi)|-,m 2 ) = - -|+,mi)|-,m 2 ) (96) 

HG\-,mi}\+,m 2 ) = - — |-,mi)|+,m 2 ) (9c) 

H G \-, mi )\-,m 2 ) = -~(h+I 2 +2)\-,m 1 )\-,m 2 ), (9d) 
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Figure 3. The different values of mi and 7712 arranged on a grid. The diagonal 
lines mark the states with constant magnetization. The dotted bended lines 
connect the states with interchanged magnetization, combined in our ansatz 
(11) . The solid bended lines connect the states which are combined in order 
to construct the set of linearly independent eigenstates spanning the eigenspace 
to the eigenvalue (J — 2A)/i. 

where we abbreviated \I ± \,mi) =: |±,mj). Obviously, the states |±, mi) |=F, m^) 
are degenerate with the eigenvalue being independent of /. As seen above, the highly 
degenerate eigenvalue in the subspace Hd is (J e x — 2A)/A. Thus, eigenstates of H 
with this eigenvalue can be constructed by simply combining the states |±, mi)|=F, m 2 ) 
to triplet states with respect to the two central spins, meaning that they lie in the 
kernel of ([5| . At this point it is of course not clear that all eigenstates with the above 
eigenvalue are resulting through this approach. However, we will see that this is indeed 
the case. 

In other words, our goal is to eliminate singlet contributions from suitable linear 
combinations of the states |±,m 1 )|=F,W2). To this end we use an ansatz already 
accounting for the conservation of J z and the parity symmetry by superimposing 
states of the form 

\+,m)\-,M-m) ± \-,M -m)\+,m) , (10) 

where M is the eigenvalue of J z . All considerations will focus on M > 0, as states 
with M < result simply by reversing every spin. In the following analysis one needs 
to distinguish the four different cases depending on whether M is even or odd, and 
/ is integer or half-integer. This case-by-case procedure can be nicely encapsulated 
and simplified as follows by introducing i = 21 — M with i = 0, . . . , 21: In Fig. [3] 
the possible values of mj and m 2 are arranged on a grid. The diagonal lines mark 
the states of constant magnetization AI = 21 — i, where we refer to the maximal 
value on such a diagonal as m max . Obviously, we have m max = 1—1/2 for i = 21 
and m max = 1+1/2 otherwise. Following a line of constant magnetization starting 
from m max , one recognizes that from a certain value m = m m i n on, all occurring 
states result from those with larger values of m by interchanging the respective 
magnetizations (m, M — m). In Fig. [3] these "complementary" states are connected 
by dotted bended lines. It is easy to see that if i is odd, we have m min = (21 — i)/2, 
whereas for an even value of i we have to add (1/2) so that m min = (21 — i + l)/2. It 
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is now a simple fact that there are states which do not have a complement. This is the 
case for the states with to = m max if i ^ 27 and for those with m — m min provided i 
is odd or equal to zero. 

With respect to later considerations it turns out to be more convenient to use 
an ansatz which is a sum over pairs of complementary states, rather than a direct 
superposition of the states (10). Hence we introduce coefficients a m ,a' m for any state 
with m > m m i n and its complement and combine them to a sum running from m m i n 

tO Tfl rr)n ,'r- 



l±,i) 



E 



|+,to)|-,27 - i- to) ± |-,27- i - m)\+,m)) 



+ 9(m)a' m (| + , 2I-i- to) |- m) ± |-, m)\+, 21 - i - to))] (11) 

The complements of the respective to = m max states automatically vanish, whereas 
the function 9(m) accounts for the m m i n states without a complement: 

8(m) = 6(to - m min - 1) + 6. mod 2fi S m , mmin (12) 

where the Heavysidc function Q(x) is unity for any x > and zero otherwise. 
Clearly, the ansatz (11) is an eigenstate of Hq, 

H a \±,i)=~\±,i), (13) 
consisting of triplet and singlet terms. Eliminating the latter by demanding 

(a m /i 4 '(to) (21 — i — to) =f Q(m)ot m n~ (to)/x~ (2/ — i — to)) 



E 

m— m„ 



to ,21 — i — m 

2' 



=F 127 — i - to - 



-, TO 

2' 2' 



+ (a m /Lt (to)/U (2/ — i — to) =p 8(m)a' m ii + (m)ii + (21 — i 



2' 



21 



T|27- 



(14a) 

to)) 

= 0,(146) 



we arrive at an eigenstate of i7 c . 

Let us first consider the two particularly simple cases i — and i = 21. For 
i = 0, i.e. M = 27, the sum consists of only one term to = 7 + 1/2. In this case 



the contributions related a' I+1 ^ 2 to in (14a) and (146) are vanishing, implying that for 
positive parity the unwanted singlet terms are automatically zero. This means that the 
largest degenerate multiplet with J = 21 always has positive parity, as demonstrated 
by our numerics. 

In the other case i = 27, i.e. M = 0, one easily sees that 

H + (m)n + (—m) = /i~ (m)fi~ (—m) . (15) 

If i is even, this condition means that for every to the singlet terms can be eliminated 
by simply choosing a m = ±a' m . Therefore, in this case we always have an equal 
number of multiplcts with positive and with negative parity. As mentioned above, for 
an odd value of i the summand with to = m m i n = does not have a complement. 
However, for positive parity the unwanted terms vanish automatically so that the 
number of positive multiplets is larger by one than the number of negative multiplets. 
In total, we get 27 solutions as suggested by our numerics. 
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These solutions result from demanding that the terms in (14a) and (14a) vanish 
separately, while, strictly speaking, only their sum is required to be zero. However, it 
is indeed simple to see that there are no further solutions: Demanding that the total 
sum vanishes leads to the conditions 

fi + (m)fi + (-m) (a m T a' m ) 

= ±/i + (m - l)^ + (-rn + 1) (a m _i T a4i-i) i 

and (^ai =p a'i J = (aj_i TQ^i^ = 0, which obviously give the same solutions as 

above. In summary, the resulting eigenstates at i = 21 (M = 0) can be formulated 
most compactly as 

|±,2/,m) :=|tt) (\m-~,-m-~)±\-m-±,m-~) 

- Ill) (\m + \,-m + ~)±\-m+^ m + . (16) 

That He lies /tiZZj/ in the kernel of ^ becomes clear at this point: There are 21 
degenerate multiplets with alternating parity, each of which gives one state with 
M = 0. Above we constructed states which are superpositions of exactly those states 
and lie in kernel of 15ft. They can be combined to give eigenstates of J 2 , so that %d 



can be constructed simply by applying J . From 



J, H c 



J, H r 



= it follows 



[J ± ,[H G ,H C }]=0, (17) 

meaning that a state resulting from the application of J to a state lying the kernel 
of (§ again lies in the kernel of Therefore the full degenerate subspace is located 
there. 

3.2. Complete construction 

Now we come to the construction of the full degenerate space Hp . In an immediate 



approach we follow the route described above and combine the states ( 16 ) to 
eigenstates of J 2 such that Wd can be generated by applying J ± . Unfortunately 
the construction of J 2 eigenstates is possible only up the solution of a homogeneous 
set of equations with a (symmetric) tridiagonal coefficient matrix, which has to be 
carried out numerically. Due to the simple shape of the matrix such a problem has 
the very low complexity of 0(2/) so that even systems of realistic size with respect 
to experimental situations in for example semiconductor quantum dots / ~ 10 6 can 
be treated on conventional computers [TJ [21 [3l [4]. However, in a second approach 
we construct a basis of Hd in a fully analytical fashion. The resulting basis states 
are eigenstates of J z and r, but they are neither orthogonal nor do they satisfy the 
J 2 symmetry. Nevertheless, for both applied as well as more mathematical future 
considerations it will be helpful to have closed analytical expressions at hand. 

3.2.1. First approach: Construction of eigenstates of J 2 with M = As mentioned 
above, ou r fi rst approach consists in using the particularly simple solutions for i — 21 



operators J 



— K 

given in ( 16 1 by combining them to eigenstates of J such that applying the ladder 
generates the full space H d ■ Hence we demand 

J 2 M±,2I,m) = J(J+1) Pm\±,2I,m) 
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^J 1 A n |±,2/,m)-J(J + l) Pm\±,2I,m): 

m—m m i n rn—iTirnin 

Explicitly this reads 

M 

Iff) E (2^+l)-2(m+-)(m--)-J(J+l) 



x Mm 



1, 



2' 
1 



-m 



± 



1 

2" 



+ z/ + (m — -)v (—to — — ) |to + 



+ z/ (to — -)j/ + (— to — — ) ( 1 771 



± 

± 



, TO - 



-, m - 



(18a) 
(186) 
(18c) 



-III) E & 



2J(/+l)-2(m+i)(m-i)- J(J+1) 



, 1 1. 1 1. 

x I \ m + ~ m + 2' ± I _ m + 2'™ + 2' 



1 



1, 



+ z/ + (mH — )v ( — to H — ) |to + 



± 



-, m 



+ v (to H — )v + (— to ) I |m , —to H — ) ± | — to H — , m ) 



0, 



where i/^ (to) = ^/ /(/ + 1) — to(to ± 1) and hence 

i/+(m- )j)v~{-m - i) = !/ + (m + i)^(-m+ i) 

v~{m- ^)v + {-m- i) = i/ _ (m+ ^)f + (-m- i), 

which is plausible, because the | ft) an d I II) terms must vanish separately. Note 
that all components with | fl), | It) are equal to zero. It is now simple to see that the 
state in ( |186| ) for some to is identical to the one in ( 186) for (to + 1) and to the one in 
(18c) for (m + 2). For an even value of i eliminating these terms gives the following 
set of equations 



v + (m )u (—to ) 

x 2' K 2' 



+ Pm+l 



+e(I - - - m)f3 m+2 



21(1 + 1) - 2(to + |)(m + \) - J( J + 1) 
v + (m + to — — ) =11 



/3i [21(1 + 1) - J(J + 1) ± /(/ + 1)] + /3 f [^-(l> + (-2)] = 0, 

where to = 1/2, — 3/2. This yields a symmetric tridiagonal matrix. However, 
the symmetry of the matrix is destroyed if i is odd. In this case we have 



(3 m (l ± S m fi) 



v + (m )u (—to ) 

v 2' K 2' 

3 

+6(1 - - - m)/3 m+2 



■A 



m+l 



2J(J + l)-2(m+|)(m+i)- J(J+1) 



^+(to+ (-to - ^) 



= 



2J(/+1)- J(J+1)- 



1 







+ ft 





= o, 
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where m = 0,...,/ — 1/2. The two above systems now have to be solved numerically 
for the different values of J. 



3.2.2. Second approach: Explicit elimination of singlet contributions via ansatz 
Our second approach, which in contrast to the above one will lead to closed 
analytical expressions for the degenerate eigenstates, consists in directly determining 



the constants a m , a' m for every value of M. As already used above, considering (14a 



for some m and (146 1 for (m— 1), one sees that the respective states become identical 
up to a factor (^1). The idea is now to eliminate these terms systematically, so that 
we get a sufficient number of linearly independent eigenvectors. As indicated in Fig. 
[3] by the solid bended lines, this can be done by simply superposing an increasing 
number of successive terms and choosing all other constants to be equal to zero. Of 
course these solutions are by no means unique. We just choose the most compact 
ones. For an odd value of i this yields the following quite cumbersome solutions 



= (-l) K - A (= F l)' i 



t-A-l 



N K 



where A = 0, . . . , (re 



2 

1) and 



' -A)/z-(J-i- 



•A) 
(19) 



N K 



(J-i-i+w)- 



(M + (/-^)) 2 Ton/^)) 2 



\- K )y+ (J-i-j + «))' 

fj,- (I + | — re)/Lt~ (/ — i — | + re) 



where the first line refers to re = 1, ...,(£— l)/2 and the second line applies to 
k = (i + l)/ 2. For even i and negative parity the solution coincides with the first 
line of (3.2.2), where now re = 1, . . . ,i/2. Considering positive parity we get: 



(i+i- K )v-(i. 



4-/r i In +/r * 1\ -/r i In -/t t 1- 



with re = 1, ...,(£ — 2)/2 for the first line and re = z/2 for the second one. Due to 
the presence of an a' m . term, in contrast to the case of an odd i, here we have an 
additional solution. This results by simply choosing all constants to be equal to zero 
except for oe mmin and a' m . , which are determined by eliminating the (14&) term: 

l > 



2 ' 2 



I) 



J Oil 



- - ^(/-| + |) M+ (J-|-|)^-^ ^ ^ (20) 
Note that all the remaining constants are determined by the normalization condition. 
Let us give a brief discussion of the above results. With respect to subspaces of fixed i, 
the degeneracies shown in Figs. [T] and [2] yield the pattern shown in Tab. [T] Obviously 
for any i there are (i + 1) states. If i is odd, there is an equal number of states with 
positive and with negative parity, whereas for an even value of i the number of states 
with positive parity is larger by one than the number of states with negative parity. 
This is perfectly reproduced by our solutions: For an odd i the index re in (3.2.2 ) runs 
up to (i + l)/2 for each parity, meaning that there are (i + 1) solutions in total. If i 
is even, (3.2.2) yields i/2 solutions for each parity and an additional one for positive 
parity. 
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Table 1. Numerically detected degeneracy pattern. 



4. The role of the inversion symmetry 

In the preceding section we have constructed the full degenerate subspace by 



determining the coefficients in (111 so that we arrive at triplet states of the two 
central spins. Obviously, such a construction is still possible if inversion symmetry 
is broken. Note that if Ji ^ I 2 , additional labels for the spin length have to be 
introduced in ([8| and (9). However, it is simple to see that in general our states are 
no longer eigenstates with respect to H G , because the degeneracy between the Hq 
eigenstates |±, mi) |=F, m 2 ) is lifted. Indeed this can be easily recovered by demanding 
E + (Ai,h) + E_(A 2 ,I 2 ) = E-(Ai,h) + E + (A 2 , h), yielding the quite remarkable 
relation 

^idim(Wi) = A 2 dim(H 2 ), (21) 

where dim(Kj) = 2/; + 1 denotes the dimension of the Hilbert space associated 
with Ij. Note that the bath spins Y^jLiH couple to different values of Jj so that, 
by a corresponding choice of Ai,A 2 , several different degenerate subspaces can be 
constructed. 



Relation (21 1 means that the inversion symmetric case is only an example of a 
whole class of systems exhibiting the same type of systematic degeneracy. In Fig. [4] 
we plot the relevant part of the spectrum for 1\ =^ 1 2 with I 2 > I\ ■ In the upper panel 



the couplings violate (21) and consequently the degeneracy between the multiplets 



is lifted. In the bottom panel it is recovered by choosing Ai,A 2 according to (21), 
leading to 

A > = jTTjT+T' 4 (22o) 

A > = ST5TT A < 22 » 

In direct analogy to the inversion symmetric case the branch begins at (ii + I 2 ) and 
ends at (I 2 — I\ + 1). 



The relation (21) has a concrete physical meaning: Consider a semiconductor 
double quantum dot. Here the electron spin interacts with the surrounding nuclear 
spins via the hyperfine interaction, yielding a system of two coupled Gaudin models. 
The role of the couplings Ai,A 2 is played by the overall coupling strengths of the 
respective dots, given the sum of all hyperfine coupling constants (which depend on 
the properties of the respective material) . The size of the Hilbert spaces results from 
the spatial extent of the respective electron wave function. If it is e.g. stretched 
over a larger area, each individual coupling decreases, but the sum remains unaltered. 



Hence, in an approximative sense, the relation (21 ) can always be realized by properly 
adjusting the electron wave function. 



Unexpected systematic degeneracy in coupled Gaudin models 

20 1 1 1 1 1 1 . 1 1 



12 



16 - 



< 12- 



12 



UJ 



I, =20; I 2 =35 
J =50A 



A =1.27 A; A,=0.73A 



21) 



41) 



51.) 



hi) 



Figure 4. Formerly degenerate branch of the spectrum for I = 20 and J = 25A. 
The inversion symmetry is broken by choosing I\ == Ii. The deviation from 
the degenerate case is stronger for multiplets with a small value of the quantum 
number J than for those with a large value. The bottom panel shows the spectrum 
with the degeneracy recovered by choosing A\ == A2 due to (22). 



With respect to possible future applications of Hd it is important to note that 



for parameters only weakly violating (21 1, the multiplets are still nearly degenerate: 
Let us fix Ai 



A 2 and vary I\,l2 so that ( |21[ ) is violated. From the eigenvalues of 
Hq it is clear that the degeneracy is lifted in a continuous way. Furthermore, as can 
be seen very well in the upper panel of Fig. [4j the influence of I\ 7^ I2 on multiplets 
with small quantum numbers J is much stronger than on those with large values of 
J. This is also the case if we analogously choose A± 7^ A 2 . 



5. Conclusion 

In summary we have reported an unexpected systematic degeneracy in an inversion 
symmetric system of two coupled Gaudin models with homogeneous couplings. This 
leads to a degenerate subspace of macroscopic size. We have constructed the complete 
degenerate subspace, which is fully located in the kernel of the commutator between 
the two Gaudin models and their coupling term. Furthermore we have studied the 
role of the inversion symmetry. Indeed it turns out that the inversion symmetric 
case is only an example for a whole family of systems all of which share the same 
type of systematic degeneracy. This exclusively originates in the degeneracy of two 
eigenspaces of the Gaudin part of the Hamiltonian, yielding a remarkable relation 
between the dimension of the bath Hilbert spaces and the couplings. 

Nevertheless, so far we have not been able to detect the (possibly continuous) 
symmetry underlying this remarkable degeneracy, i.e. a set of generating operators 
that would connect the highly degenerate multiplets. This question remains as an 
important but probably rather intricate problem for further studies. Furthermore it 
would be fruitful to study applications of the degenerate space especially in the context 
of solid state quantum information processing. 
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